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This  report  examines  details  of  a  one-dimensional  (ID)  atmospheric  bound¬ 
ary  layer  model  to  establish  the  proper  functioning  of  its  soil,  plant,  and  at¬ 
mospheric  physics.  To  achieve  this  goal,  I  inspect,  repair,  and  modify  a  com¬ 
puter  program  that  scientists  at  the  Hebrew  University,  Department  of  Soil 
and  Water  Sciences,  gave  to  me  years  ago.  The  computer  program  was  exer¬ 
cised  to  determine  if  the  model  results  are  stable  when  initial  conditions  are 
changed  and  to  determine  whether  the  results  are  sensible  and  generally 
consistent  with  observed  data.  To  show  this,  I  present  a  time  series  of  the 
modeled  surface  energy  budget  and  modeled  profiles  of  boundary  layer 
wind  speed,  potential  temperature,  and  specific  humidity  for  daytime  (at¬ 
mospherically  unstable  conditions)  and  for  nighttime  (atmospherically  sta¬ 
ble  conditions).  I  compare  these  results,  wherever  practical,  with  observed 
meteorological  data.  From  these  results,  I  infer  how  well  the  transfers  of 
momentum,  heat,  and  moisture  from  one  model  layer  to  the  next  are  char¬ 
acterized.  I  also  present  root  mean  square  error  and  d  values,  where  d  is  an 
index  of  agreement,  to  summarize  the  model  results  and  comparison  with 
observed  data.  From  the  results,  I  find  that  the  ID  model  is  functioning 
properly  in  solving  for  many  parameter  relationships  and  is  as  reliable  as 
the  earlier  models  of  this  type  in  predicting  the  general  features  of  bound¬ 
ary  layer  development. 
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1.  Introduction 


The  atmospheric  boundary  layer  is  also  called  the  friction  layer.  It  extends 
from  the  earth's  surface  to  the  geostrophic  wind  (or  gradient  wind)  level. 
Above  this  layer  is  the  free  atmosphere — ^where  the  frictional  influence  of 
the  earth's  surface  is  greatly  diminished,  allowing  for  an  approximate  bal¬ 
ance  between  pressure  gradient  and  Coriolis  accelerations  (Huschke,  1959). 
The  daytime  boundary  layer  is  often  observed  to  heights  of  1  to  3  km 
above-ground  level  (agl),  generating  convective  eddies  (thermal  updrafts 
and  downdrafts)  and  relatively  well-mixed  profiles  of  wind  speed,  tem¬ 
perature,  and  moisture.  In  contrast,  the  nighttime  boundary  layer  is  char¬ 
acterized  by  a  temperature  inversion  caused  by  a  strong  radiative  cooling 
at  the  surface.  The  nighttime  boundary  layer  can  be  at  heights  of  300  to 
400  m  agl.  A  low-level  wind  maximum  or  jet  can  sometimes  develop  in  the 
nighttime  boundary  layer  with  faster,  warmer  air  flow  aloft,  and  slower, 
cooler  air  at  the  surface.  This  development  of  a  low-level  jet  can  promote 
even  further  cooling  of  the  surface  layer  air,  unless  the  inversion  breaks 
down  or  overturns,  which  may  occur  because  of  increased  wind  shears 
or  other  larger-scale  instabilities,  perhaps  also  in  combination  with  atmos¬ 
pheric  waves  (Businger,  1973). 

In  formulating  an  atmospheric  boundary  layer  study,  one  needs  a  useful 
model  calculation.  That  is,  are  the  soil,  plant,  and  atmospheric  physics 
in  the  model  functioning  correctly  and  reliably?  To  achieve  this  goal,  one 
needs  to  study  the  details  of  a  model  program  to  determine  if  the  model 
results  are  stable  and  whether  the  results  are  sensible,  i.e.,  generally  consis¬ 
tent  in  comparison  with  observed  data. 

In  this  report,  I  use  a  computer  model  (see  sect.  2)  that  scientists  at  the  He¬ 
brew  University,  Department  of  Soil  and  Water  Sciences,  gave  to  me  sev¬ 
eral  years  ago.  To  test  the  model,  I  generate  a  time  series  of  the  surface 
energy  budget  and  profiles  of  boundary  layer  wind  speed,  potential  tem¬ 
perature,  and  specific  humidity  for  daytime  (atmospherically  unstable  con¬ 
ditions)  and  for  nighttime  (atmospherically  stable  conditions).  I  compare 
these  model  results,  wherever  practical,  with  observations  from  days  33 
and  34  of  the  Wangara  experiment  (Clarke  et  al,  1971).  From  these  results, 
I  examine  how  the  model  characterizes  transfers  of  momentum,  heat,  and 
moisture.  I  attempt  to  quantify  these  results  by  preparing  root  mean  square 
error  (rmse)  and  d  values,  where  d  is  an  index  of  agreement  (described  in 
sect.  5),  to  summarize  the  model  comparison  with  observed  data. 
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2.  Model  Description 


The  one-dimensional  soil,  plant,  and  atmospheric  model  used  in  this  report 
has  been  documented  previously  in  Pielke  and  Mahrer  (1975),  McNider 
and  Pielke  (1981),  and  Avissar  and  Mahrer  (1988).  It  is  a  first-order  clo¬ 
sure  model  to  calculate  the  transfer  of  momentum,  heat,  and  moisture  at 
the  surface  and  aloft.  It  uses  an  implicit  finite  difference  scheme  to  inte¬ 
grate  the  boundary  layer  and  soil  diffusion  equations.  It  contains  a  com¬ 
plete  model  of  the  surface  energy  budget  and  a  time-dependent  calculation 
of  the  height  of  the  daytime  planetary  boundary  layer.  I  added  the  formu¬ 
lation  suggested  by  Smeda  (1979)  for  the  time-dependent  calculation  of  the 
height  of  the  nighttime  stable  layer.  The  model  surface  layer  turbulence 
scaling  is  as  described  by  Zilitinkevich  (1970)  and  Businger  et  al  (1971). 

As  initial  input,  the  model  requires  day  of  the  year  and  time  of  day  data, 
latitude  and  longitude,  fraction  of  sky  cloudiness,  and  ground-cover  data 
(i.e.,  canopy  height,  leaf  area  index,  surface  reflectivity  (albedo),  thermal 
emissivity,  and  surface  roughness).  Typical  values  for  albedo  and  surface 
roughness  have  been  documented  by  Hansen  (1993a,  1993b).  The  model 
computes  soil  properties  (e.g.,  hydraulic  and  thermal  conductivity  and  soil- 
specific  heat  capacity)  from  inputs  of  soil  water  content,  porosity,  texture 
(i.e.,  proportion  of  sand,  clay,  and  organic  matter),  and  bulk  density.  Other 
model  constants  include  subsoil  properties,  such  as  plant  root  density  and 
distribution.  Also,  the  model  requires  an  initial  profile  of  wind  speed,  air 
temperature,  and  specific  humidity  (atmospheric  water  vapor  content)  from 
the  ground  level  to  the  top  of  the  model.  Model  parameters  for  this  report 
are  summarized  in  table  1. 


The  model  equations  for  the  time-dependent  calculation  of  the  winds  (u 
and  V  components),  potential  temperature  {0),  and  specific  humidity  {q) 
over  flat  earth  can  be  expressed  as 


du 

dv 

de 

~dt 

dq 

~di 


/(%-h)  +  T(a„A), 


(1) 

(2) 

(3) 

(4) 


where  /  denotes  the  Coriolis  parameter  (the  deflecting  force  caused  by 
the  earth's  rotation  acting  upon  moving  air),  the  subscript  g  refers  to  the 
geostrophic  wind  (such  that  the  first  term  on  the  right  contains  the  pres¬ 
sure  gradient  acceleration).  Km  denotes  the  eddy  transfer  coefficient  for 
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Table  1 .  Boundary  layer 
model  parameters. 


K,  =  K,,  + 


Parameter 

Value 

Number  of  vertical  levels 

32  (2, 10, 50, 100, 150, 200,  250, 300, 
350, 400, 450, 500,  550,  600,  650, 

700,  750,  800,  850,  900,  950, 1000, 
1100, 1200, 1300, 1400, 1500, 

1600, 1700,  1800, 1900,  2000  m) 

Latitude,  longitude 

34.50  S,  144.93  E 

Surface  roughness 

0.0045  m 

Vegetation  (sparse) 

0.01m 

Surface  albedo 

0.20 

Surface  emissivity 

0.98 

Soil  water  content 

0.08  m®/m® 

Average  soil  density 

1600  kg/ m® 

Soil  texture 

(%  sand,  %  clay,  %  organic) 

28.0,  70.0,  2.0 

Day,  month,  year 

16-17  August  1967 

Geostrophic  wind  (ug,  Vg) 

—5.34  m/ s,  —0.43  m/s 

Initial  time 

09  It 

Time  step 

10  s 

momentum,  and  Kh  =  Kq  denotes  the  transfer  coefficients  for  heat  and 
moisture.  In  the  surface  layer,  and  are  calculated  as  Km  =  ^u^z/cpm 
and  Kh  =  \M^z/4>h,  respectively,  where  k  is  Karman's  constant  (=  0.4);  z  is 
height  above  ground  level  (in  meters);  u*  is  the  friction  velocity  (in  units  of 
m/s“^),  which  relates  to  surface  stress;  and  the  ^m  and  (ph  are  nondimen- 
sional  lapse-rate  functions,  which  account  for  surface-layer  stabilities  other 
than  neutral.  A  list  of  symbol  definitions  is  provided  in  the  appendix. 

In  the  unstable  daytime  boundary  layer,  the  model  derives  the  transfer  co¬ 
efficients  as  a  function  of  height,  as  suggested  by  O'Brien  (1970).  This  func¬ 
tion  can  be  expressed  as 


(  {z,-zf\ 


K.,  + 


Z^)  X 


dz 


-b 


{Zi  -  Zs) 


(5) 


where  z  is  height  above  ground  level  and  Zi  and  Zg  refer  to  the  heights  of 
the  top  of  the  surface  layer  and  top  of  the  boundary  layer,  respectively.  The 
profile  function  is  applied  similarly  for  the  momentum,  heat,  and  moisture 
coefficients.  The  model  sets  Zg  =  O.OTz*,  Kz^  =  1.0,  and  Kz  =  1.0  for  z  >  Zi. 
Figure  1  illustrates  the  geometry  (i.e.,  vertical  levels)  of  the  model. 
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Figure  1.  Sketch  of 
geometry  (i.e.,  vertical 
levels)  of  model. 


i  =  N  (model  top) 
j  =  N-^ 
j  =  N-2 


Kfj,  Ri 
u,  V,  6,  q 

Boundary  layer 


z  =  0^ 


Soil  layer 


7  +  1 
.7  +  1/2 


Soll-plant-atmosphere  Interface 
.7-1/2 

^  7  +  1/2 

i  =  N  (model  bottom) 


The  growth  of  the  daytime  planetary  boundary  layer  is  calculated  with  an 
expression  derived  by  Deardorff  (1974),  which  can  be  written  as 


dzj 

dt 


-  W,i  = 


1.8  {w*  +  l.lu?  —  S.SuifZi) 
9{i)  {^)+9w^*+7.2u^ 


(6) 


where  Wzi  is  the  vertical  velocity  at  z*  (wzi  is  assumed  to  be  negligi¬ 
ble  or  equal  to  zero,  since  time-dependent  calculations  of  vertical  veloc¬ 
ities  are  not  explicitly  derived),  dO'^  jdz  is  the  vertical  gradient  of  poten¬ 
tial  temperature  in  the  stable  air  immediately  above  the  boundary-layer 
top,  6s  is  the  potential  temperature  at  the  top  of  the  surface  layer  Zs,  and 

=  [{—g/9)  u^O^Zi]  ^3  is  the  vertical  velocity  scaling  variable,  an  implicit 
calculation  of  buoyancy,  where  g  is  acceleration  caused  by  gravity,  and 

~  tfif'  1-^®  potential  temperature  scaling  constant. 

During  the  nighttime,  when  the  surface  layer  is  stable  (i.e.,  0*  >  0),  the 
model  derives  eddy  transfer  coefficients  above  the  surface  layer,*  as  sug¬ 
gested  by  Blackadar  (1979),  i.e.. 


Km{z)  =  Kh{z) 


sP{l  -  18Ri)°•^Ri  <  0 
^  si  (Ricrit  Ri)  /RicritjO  <  Ri  ^  Ricrit  ’ 
0,  Ri  >  Ricrit 


(7) 


where  s  is  the  local  wind  shear,  s  =  y  {du/dz)"^  -F  {dv/dz)‘^  ,  and  Ri  is  the 
ratio  of  thermal  to  mechanical  (wind  shear)  production  turbulent  energy 
called  the  Richardson  number,  so  that  Ri  =  f  §f/  [{du/dz)^  -F  {dv/dz)^^. 

‘Note:  Za  =  0.15  Zi  while  0*  >  0. 


4 


Ricrit  is  the  limiting  value  of  the  Richardson  number,  and  it  is  often  as¬ 
sumed  that  Ricrit  =  0.25,  even  though  Ricrit  =  1-0  can  sometimes  be  a  use¬ 
ful  approximation  (Avissar  and  Mahrer,  1988).  The  length  I  (in  meters)  is 
generally  thought  of  as  the  width  of  turbulence  and  can  be  characterized 
by  the  formulation  reported  in  Blackadar  (1979)  for  z  >  Zg  as 

^  +  o.ooesrt*//)  ■ 

The  height  of  the  nighttime  planetary  boundary  layer  is  calculated  with  the 
formulation  suggested  by  Smeda  (1979),  which  can  be  expressed  as 


dzi 

dt 


11^ 

0.06^ 

ZiJ 


1  - 


S.Szif 


(9) 


The  surface  energy  budget  is  applied  to  the  soil,  the  plant  canopy  and  the 
air  throughout  the  canopy,  and  the  thin  layer  of  air  that  extends  above  the 
canopy  top  (Avissar  et  al,  1986).  Air  temperature  and  water  vapor  content 
for  each  layer  are  derived  because  of  energy  and  water  vapor  transfers 
within  the  system.  The  net  radiative  flux  is  approximated  as  functions  of 
transmission,  albedo,  leaf  area,  and  soil  wetness.  The  model  radiation  and 
energy  budget  as  described  in  Pielke  (1984)  are 

{1  -  A) Rsi  +  Rli  -  Rli  -  pCpU^O^  -  pL^u^q^  +  Qs  =  F,  (10) 

where  A  is  surface  reflectivity  (albedo);  Rg^,  R^i,  and  Rlj  are  the  incom¬ 
ing  solar,  incoming  long-wave,  and  outgoing  long-wave  radiative  fluxes, 
respectively;  p  is  air  density;  Cp  is  the  specific  heat  of  air  at  constant  pres¬ 
sure;  u^,  9^,  and  g*  are  the  surface-layer  turbulence  scaling  parameters  for 
wind  speed,  temperature,  and  moisture,  in  that  order;  is  the  heat  of 
transformation  for  water  vapor;  Qg  is  the  soil  heat  flux;  and  F  is  the  func¬ 
tion  applied  in  solving  for  the  surface  temperature  Ogic-  Using  the  Newton- 
Raphson  iterative  algorithm,  I  can  approximate  each  successive  estimate  of 
the  surface  temperature  as 

(11) 


The  model  equations  for  the  soil  layer  and  plant  canopy  are  described  in 
much  greater  detail  by  Avissar  and  Mahrer  (1988),  with  one  exception.  The 
model  formulation  for  downward  long-wave  radiation  is  instead  based  on 
the  empirical  relationship  reported  by  Paltridge  and  Platt  (1976),  which  can 
be  expressed  as 


Rli  =  -170.9  +  lA95aT^ +  0.3c\decaT^,  (12) 

where  a  =  5.6697  x  10“®  W/m^  deg“^,  Tr  is  the  reference  level  (~2  m) 
temperature  in  Kelvin,  cld  is  the  cloud  amount,  Sc  is  the  emissivity  of  the 
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cloud  base,  and  Tc  is  the  temperature  of  the  cloud  base  in  Kelvin.  This  for¬ 
mulation  is  much  simpler  to  apply  than  computing  the  upward  and  down¬ 
ward  long-wave  radiation  according  to  the  concentrations,  path  lengths, 
and  emissivities  for  water  vapor  and  carbon  dioxide.  Also,  this  expression 
was  used  successfully  in  two  previous  studies  (Rachele  and  Tunick,  1994; 
Tunick  et  al,  1994). 
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3.  Numerical  Methods  and  Boundary  Conditions 


The  numerical  integration  of  equations  (1)  to  (4)  is  achieved  with  the  use  of 
a  generalized  form  of  the  Crank-Nicholson  implicit  finite  difference  scheme 
(Paegle  et  al,  1976)  to  solve  for  the  vertical  transfer  of  momentum,  heat,  and 
moisture  in  the  atmosphere.  The  diffusion  terms  in  equations  (1)  to  (4)  can 
be  expressed  in  terms  of  both  a  current  r  and  future  r  +  1  timestep: 


^r+l  _ 


At 


-^i+1/2  ^^j+i/2 

1  1/2 


(13) 


where  (5r  +  (5t+i  =  1,  Azj  =  Azj+i  =  Zj+i  -  Zj,  and 

Azj-i  =  Zj  —  Zj-i.  In  this  study,  j3r  =  0.75  and  Pt+i  =  0.25  are  the  weights 
to  the  current  and  future  contributions  to  the  numerical  approximation 
(Pielke,  1984).  Equation  (13)  can  be  rewritten  as  aj4>j-i  +  bj4>j  +  Cj4>j+i  =  dj 
such  that 


PT+iKj_i/2At 

AzjAzj_ii2 


—  (Pn  + 


l3rKjj^ll2At  ^ 


AzjAz 


^  ^  PT+lKj_^_lj2^t  ^  I^T+lKj_l/2At 
AzjAzj+ij2  AzjAzj_ii2 

0J+1  -  (5rKj_y2^t 

i+1/2  ^ZjAZj_i/2 


T+1 

j 


PT+iKj_^_ij2At 

AzjAzj+i/2 

(14) 


The  model  solves  equation  (14)  using  the  coefficients  aj,  bj,  Cj,  and  dj  and 
the  algorithm  described  by  Ahlberg  et  al  (1967),  which  can  be  expressed  as 

(t)j  =  Xj  +  yj(j)j+i,  (15) 

where  Xj  =  {dj  —  ajXj-i)/pj,  pj  =  bj  +  ajPj-i,  and  pj  =  —Cj/pj,  and  where 
xo  =  yo=  0.  When  the  pressure  gradient  and  Coriolis  accelerations  in  equa¬ 
tions  (1)  and  (2)  are  included,  the  coefficient  dj  (i.e.,  the  right-hand  side  of 
eq.  (14))  becomes  dj  =  dj  +  f  {v  —  Vg)  At  and  dj  =  dj  +  f  {ug  —  u)  At, 
respectively. 

The  boundary  condition  applied  at  the  top  of  the  model,  i.e.,  j  =  N,  for  this 
solution  is  (j)N  =  0.754>n-i  +  0.25(I)n-2,  so  that  the  variables  near  the  top  of 
the  model  are  not  fixed  (i.e.,  not  held  to  their  initial  values).  This  condition 
appears  to  behave  well. 
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4.  Days  33  and  34  Observations 


The  Wangara  experiment  data  were  collected  over  a  large  flat  area 
60  X  60  km),  consisting  mainly  of  thin  low  bushes  and  desert  grasses. 
The  field  study  was  conducted  near  the  town  of  Hay,  New  South  Wales, 
Australia  (34.50  S,  144.93  E),  during  the  1967  southern  hemisphere  winter 
(Clarke  et  al,  1971).  Five  types  of  data  were  collected:  (1)  hourly  meteoro¬ 
logical  tower  measurements  of  wind  speed,  specific  humidity,  and  temper¬ 
ature;  (2)  hourly  surface  measurements  of  the  net  radiative  and  soil  heat 
flux;  (3)  measurements  of  wind  speed  and  wind  direction  from  0  to  2  km 
agl  obtained  by  hourly  pibal  (pilot  balloon)  flights;  (4)  measurements  of 
pressure,  temperature,  and  mixing  ratio  (weight  of  water  vapor /weight  of 
dry  air)  to  a  height  of  2  km,  taken  by  radiosonde  (instrumented  balloon 
sonde)  flights  at  3-hour  intervals;  and  (5)  reports  of  fractional,  low,  high, 
and  total  cloud  cover  (including  cloud-type  descriptors). 

These  data  are  sufficent  to  test  the  functioning  of  boundary  layer  mod¬ 
els  (e.g.,  Pielke  and  Mahrer,  1975;  Yamada  and  Mellor,  1975;  McNider  and 
Pielke,  1981).  Of  the  1050  published  reports  (profiles),  I  used  those  data  for 
the  27-hour  period  from  09  It  (local  time)  on  16  August  1967  to  12  It  on 
17  August  1967.  This  time  period  was  characterized  mainly  by  clear  skies, 
relatively  light  daytime  wind  speeds,  dry  soil,  and  strong  surface-based 
convective  heating.  In  contrast,  the  nighttime  boundary  layer  winds  over 
this  time  period  were  quite  strong,  and  the  surface  layer  cooled  significantly. 
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5.  Model  Results 


Figure  2  gives  the  time  series  of  the  modeled  u  and  v  components  of  wind 
speed  and  the  time  series  of  the  total  wind,  at  2  m  agl,  in  com¬ 

parison  with  days  33  and  34  observations.  The  observed  wind  speeds  are 
shown  as  mostly  overpredicted,  especially  through  the  period  12  to  18 
It.  As  a  result,  the  values  of  the  index  of  agreement*  calculated  for  these 
data  are  about  50  to  60  percent  (see  table  2).  These  comparisons  might 
have  been  improved  if  adjustments  to  the  geostrophic  wind  and  thermal 
wind  (i.e.,  Ug,  Vg,  dugjdz,  and  dvgjdz)  were  included  in  the  calculation. 
I  use  zq  =  0.45  cm  based  on  modeling  the  surface  energy  budget  over  a 
barren  field  (Tunick  et  al,  1994).  Several  others,  however,  have  used  val¬ 
ues  of  zo  other  than  0.45  cm;  for  example,  Pielke  and  Mahrer  (1975)  used 
Zq  =  0.12  cm,  as  reported  on  p.  21  of  Clarke  et  al  (1971),  and  McNider  and 
Pielke  (1981)  used  zq  =  1.0  cm,  based  on  the  recommendations  of  Yamada 
and  Mellor  (1975)  and  others  that  0.12  cm  was  too  small  for  use  with  the 
surface-layer  flux  formulations  in  their  models. 


Figure  2.  Modeled  u  and 
V  components  of  wind 


speed  and  total  wind, 
a/u^  -f  at  2  m  agl. 


compared  with  days  33 


and  34  observations. 


T3 

OA 

(D 

CL  05 


Local  time 


*The  index  of  agreement  d,  as  suggested  by  Willmott  (1981),  is  calculated  as 
d  =  1-  E  (Mi-OdVE  [\{Mi-Oi)\  +  \{Oi-Oi)\]\ 

i=l  i=l 

where  Mi  are  the  modeled  data,  Oi  are  the  observed  data,  and  the  overbar  corresponds  to 
the  mean.  The  nature  of  the  index  of  agreement  is  a  simple  scale  0  to  100  percent,  such  that 
d  ^  1  as  the  model  predictions  improve.  In  contrast,  the  rmse  is  calculated  in  the  same 

N 

units  as  the  variable,  i.e.,  rmse  =  E  ~  Oi)^ 

V  i=l 
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Table  2.  Summary  of 
model  results  in 
comparison  with 
observed  data. 


Model  parameter 

rmse 

d 

n 

Total  wind  speed  @  2  m  agl,  m/s 

1.749 

0.592 

28 

u  =  wind  speed  @  2  m  agl,  m/s 

1.187 

0.616 

28 

V  =  wind  speed  @  2  m  agl,  m/ s 

2.148 

0.480 

28 

Temperature  @  2  m  agl,  K 

0.627 

0.966 

28 

Net  radiative  flux,  W/m^ 

25.098 

0.995 

27 

Soil  heat  flux,  W/m^ 

30.273 

0.955 

27 

Boundary  layer  height,  m 

218.495 

0.888 

7 

Wind  speed  profiles,  m/s 
(12, 15, 18,  21  It) 

2.838 

0.601 

120 

Wind  speed  profiles,  m/s 
(24,  03,  06,  09  It) 

5.010 

0.372 

120 

Potential  temperature  profiles,  K 
(12, 15, 18,  21  It) 

0.435 

0.992 

120 

Specific  humidity  profiles,  g/kg 
(12, 15, 18,  21  It) 

0.209 

0.988 

120 

Figure  3  gives  the  time  series  of  the  surface  and  2-m  air  temperatures.  The 
modeled  surface  temperatures  9sic  are  shown  in  comparison  with  results 
from  an  earlier  modeling  study  reported  by  McNider  and  Pielke  (1981).  The 
modeled  air  temperatures  at  2  m,  62m,  are  shown  in  comparison  with  the 
days  33  and  34  observations.  The  results  are  in  very  good  agreement  with 
these  data,  especially  for  the  daytime  hours.  They  imply  that  the  model's 
surface  energy  budget  calculations  (shown  in  fig.  4)  are  functioning  prop¬ 
erly.  Index  of  agreement  d  and  rmse  values  for  these  data  as  well  as  the  sur¬ 
face  energy  budget  are  given  in  table  2.  The  discrepancies  between  modeled 
and  observed  data  at  nighttime,  however,  may  be  related  to  initial  values 
of  surface  roughness  and  to  certain  initial  soil  properties,  such  as  soil  water 
content. 

Figure  5  gives  the  modeled  height  of  the  planetary  boundary  layer.  The 
model  formulations  for  Zi,  which  are  based  on  calculations  of  the  surface- 
layer  turbulence  scaling  parameters,  appear  to  be  in  fairly  good  agreement, 
i.e.,  d  ^90  percent,  with  the  observations,  except  perhaps  during  the  hours 
immediately  following  sunset.  The  observations  are  few  because  they  were 
estimated  from  days  33  and  34  profiles  of  6  and  q,  taken  at  3-hour  intervals 
(see  Melgarejo  and  Deardorff,  1974,  and  Pielke  and  Mahrer,  1975). 

Figure  6  gives  the  modeled  boundary  layer  wind  speed  profiles  at  3-hour 
intervals  in  comparison  with  days  33  and  34  observations.  During  the  day¬ 
time,  i.e.,  12  to  18  It,  the  profile  average  wind  speeds  (both  modeled  and 
observed)  are  shown  to  increase  through  the  layer  2;  <  1300  m.  Similarly, 
the  modeled  and  observed  profiles  of  wind  speed  at  night,  i.e.,  24  to  06  It, 
show  the  development  of  a  low-level  jet  in  the  layer  100  <  2:  <  400  m.  How¬ 
ever,  agreement  between  the  modeled  and  observed  profiles  for  these  data 
is  generally  poor  (i.e.,  d  ^  60%,  daytime,  and  d  ^  37%,  nighttime),  particu¬ 
larly  through  the  upper  layers  and  near  the  model's  top.  From  this  result,  I 
suspect  that  either  the  upper  boundary  condition  for  equations  (1)  and  (2) 
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Figure  3.  Modeled 
surface  (6»sfc) 
temperatures  compared 
with  results  from  an 
earlier  modeling  study 
reported  by  McNider 
and  Pielke  (1981)  and 
modeled  2-m  (d2m)  air 
temperatures  compared 
with  days  33  and  34 
observations. 

Figure  4.  Modeled 
surface  energy  budget 
compared  with  days  33 
and  34  observations  of 
net  radiative  and  soil 
heat  flux. 


Figure  5.  Modeled 
height  of  planetary 
boundary  layer 
compared  with  days  33 
and  34  observations. 


Local  time 

does  not  apply  and/ or  the  observed  profile  data  reflect  changes  (over  time) 
in  the  geostrophic  wind  (see  table  3),  an  effect  that  was  not  included  in  the 
calculation.  This  can  be  addressed  in  a  separate  study. 

Figure  7  gives  the  modeled  potential  temperature  and  specific  humidity 
profiles  at  3-hour  intervals  in  comparison  with  days  33  and  34  observa¬ 
tions.  The  potential  temperature  profiles  are  well  mixed  vertically  through 
the  convective  layer,  z  <  1300  m,  increasing  (on  the  average)  over  the  day¬ 
time  hours.  The  specific  humidity  profiles  are  also  well  mixed  through 
the  convective  layer  and  are  decreasing  (on  the  average)  over  the  daytime 
hours.  The  index  of  agreement  between  these  modeled  and  observed  data 
is  d  ^  0.99.  In  part,  this  high  value  of  agreement  is  because  there  is  little  or 
no  advection  of  heat  or  moisture  over  time  in  the  upper  model  layers.  For 
layer  z  <  Zi,\  can  infer  that  the  eddy  transfer  (eddy  diffusion)  coefficients 
used  in  the  model  are  formulated  reasonably  well.  Figure  8  summarizes 
these  comparisons  of  the  modeled  results  to  observed  data. 
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Figure  6.  Model  boundary  layer  wind  speed  profiles  at  3-hour  intervals  compared  with  days  33  and  34 
observations. 


Table  3.  Days  33  and  34  Total  wind  Direction 

geostrophic  wind  data.  Local  time  (m/s)  (m/s)  (m/s)  (°) 


09 

-5.34 

-0.43 

5.36 

85.4 

12 

-5.56 

1.27 

5.70 

102.9 

15 

-6.20 

0.98 

6.28 

99.0 

18 

-6.42 

-0.72 

6.46 

83.6 

21 

-5.99 

-1.93 

6.29 

72.1 

24 

-6.93 

-2.86 

7.50 

67.6 

03  (27) 

-8.02 

-3.13 

8.61 

68.9 

06  (30) 

-7.32 

-4.97 

8.85 

55.8 

09  (33) 

-7.60 

-4.72 

8.95 

58.2 

12  (36) 

-7.95 

-3.43 

8.66 

66.7 

12 


Height  agl  (m)  Height  agl  (m) 


Potential  temperature  (K)  Potential  temperature  (K) 


Figure  7.  Modeled  potential  temperature  and  specific  humidity  profiles  at  3-hour  intervals  compared  with 
days  33  and  34  observations. 
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Observations  Observations 


6.  Conclusion 


For  this  report,  I  studied  and  modified  an  atmospheric  boundary  layer 
computer  program.  I  ran  the  computer  program,  making  changes  to  var¬ 
ious  sets  of  initial  conditions.  From  this  effort,  I  found  that  the  one¬ 
dimensional  model  functions  correctly  and  appears  to  be  as  reliable  as 
any  of  the  earlier  models  of  this  type  (e.g.,  Pielke  and  Mahrer  (1975)  and 
McNider  and  Pielke  (1981))  in  predicting  the  general  features  of  bound¬ 
ary  layer  development.  I  hope  to  use  this  work  toward  conducting  an 
additional  study  on  boundary  layer  wind  shears  under  nighttime  stable 
conditions. 

I  found  this  work  to  be  very  challenging.  In  particular,  rendering  the 
Ahlberg  et  al  (1967)  algorithm  from  the  computer  code  provided  a  good 
introduction  to  numerical  methods  used  in  computer  models. 
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Appendix — Symbols  and  Definitions 


af  coefficient  of  finite  difference  scheme 
A:  surface  reflectivity  (albedo) 
bji  coefficient  of  finite  difference  scheme 
cld:  cloud  amount  (in  tenths) 

Cji  coefficient  of  finite  difference  scheme 
Cpi  specific  heat  of  air  at  constant  pressure 
d:  index  of  agreement 
df  coefficient  of  finite  difference  scheme 
/:  Coriolis  parameter 

F:  function  applied  in  solving  for  surface  temperature  through  surface 
energy  budget 

g:  acceleration  caused  by  gravity 
j :  vertical  level  of  model 
k:  Karman's  constant  (=  0.4) 

Kh'-  eddy  transfer  coefficient  for  heat 
Km-  eddy  transfer  coefficient  for  momentum 
Kq.  eddy  transfer  coefficient  for  moisture 
1:  mixing  length 

Lt,:  heat  of  transformation  for  water  vapor 
m:  total  number  of  iteration  steps  to  solve  for  equation  (10) 

Mi',  modeled  data 

n:  total  number  of  model  or  observed  data  points 
N :  total  number  of  vertical  grid  points 
Oi\  observed  data 

pf  coefficient  in  algorithm  to  solve  for  equation  (14) 
q\  specific  humidity 

g* :  surface-layer  turbulence  scaling  parameter  for  moisture 
Qs'.  soil  heat  flux 

Ri:  ratio  of  thermal  to  mechanical  (wind  shear)  production  turbulent  en¬ 
ergy  called  Richardson  number 
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Ricrit^  limiting  value  of  Richardson  number 
Rsi :  incoming  solar  radiative  flux 
Rli-  incoming  long-wave  radiative  flux 
Rl^  :  outgoing  long-wave  radiative  flux 
s:  local  wind  shear,  s  = 
t:  time 

Tc'.  temperature  of  cloud  base  (in  Kelvin) 

Tr'.  reference  level  (~2  m)  temperature  (in  Kelvin) 
w.  east-west  component  of  horizontal  wind  speed 
Ug\  east- west  component  of  geostrophic  wind  speed 
ii* :  surface  friction  velocity 
V.  north-south  component  of  horizontal  wind  speed 
Vg-.  north-south  component  of  geostrophic  wind  speed 
Wzi'-  vertical  velocity  at  Zi 

:  vertical  velocity  scaling  variable 
Xj\  coefficients  in  algorithm  in  equation  (15) 
yf.  coefficients  in  algorithm  in  equation  (15) 
z:  height  above  ground  level  (agl) 

Zi\  height  of  top  of  planetary  boundary  layer 
Zs-  height  of  top  of  surface  layer 
[3r-  weighting  function  related  to  current  timestep 
I3r+i'-  weighting  function  related  to  future  timestep 
Ec-  emissivity  of  cloud  base 
4>h-  nondimensional  temperature  lapse  rate 
(pm'-  nondimensional  wind  shear 
cpN'-  any  profile  variables  (i.e.,  u,  v,  9,  or  q) 

6:  potential  temperature 
0* :  potential  temperature  scaling  constant 
6*sfc:  surface  potential  temperature 
92m-  potential  temperature  at  2  m 
p:  air  density 

a:  Stefan-Boltzmarm  constant 
r:  indicates  current  timestep 
T  +  1:  indicates  future  timestep 
any:  overbar  denotes  the  mean 
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